General Connectivity Distribution Functions for Growing Networks 
with Preferential Attachment of Fractional Power 
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We study the general connectivity distribution functions for growing networks with preferential 
attachment (PA) of fractional power, Ili oc k°', using Simon's method. We first show that the heart 
of the previously known methods of the rate equations for the connectivity distribution functions 
is nothing but Simon's method for word problem. Secondly, we show that for the case of fractional 
a the Z-transformation of the rate equation provides a fractional differential equation of new type, 
which coincides with that for PA with linear power, when a = 1. We show that to solve such a 
fractional differential equation we need define a transidental function T(a, 6, c; z) that we call upsilon 
function. Most of all the previously known results are obtained consistently in the frame work of a 
unified theory. 



I. INTRODUCTION 

Random network models have been a prototype model 
for a complex network system for a long time 1]. How- 
ever, at nearly the end of 1990's the scale-free networks 
were discovered from studying the growth of the inter- 
net geometry and topology. After the discovery, scien- 
tists have found that many real systems such as internet 
topology, human sexual relationship, scientific collabora- 
tion, economical network, and so on show the scale-free 
network topology. [1, 0, 0, H, H, 0, Q • The most charac- 
teristic property of the scale-free network is the scaling 
property that the distribution function P{k) of numbers 
of links with degree k is given by P{k) oc fc~^ where 
2 < 7 < oo. To describe the scaling law in the networks 
there have appeared many studies, using various statisti- 
cal and mathematical methods 



In those studies Barabasi and Albert @ first introduced 
a growing network model (called the B A- model). They 
introduced the concept of growing networks with pref- 
erential attachment (PA) that is proportional to link 
degree k and the system exhibits a scaling property of 
P{k) oc k~^. In terms of time development of the con- 
nectivity distribution function, Krapivsky, Redncr and 
Leyvraz[llj and Dorogovtsev and MendesJ^] have gen- 
eralized the PA to include the more general PA that is 
proportional to fc" and to that is proportional to attrac- 
tiveness A in addition to linear PA such a,s A -\- k, re- 
spectively. Liu, Lai and Ye[l^ have given a most direct 
generalization of the BA-model to include the fc"-PA in 
the spirit of Krapivsky, Redner and Levvrazfllj. 

On the other hand, to show the scale- free nature and to 
obtain the scaling exponent of the distribution function, 
Bornholdt and Ebel[16j have presented another method 
based on the famous Simon's method [28|. They obtained 
the exponent 7 = 1 -|- 1/(1 - /3), where /3 (0 < /3 < 1) 
is the probability that a newly attaching node is exerted 
into the system, while 1 — /3 is the probability that a new 
link is attached between the existing nodes. Therefore, 



if /3 = 1 then the system must coincide with the one 
studied by Krapivsky, Redner and Leyvraz[ll]. However, 
it does not do so, because in this limit the former provides 
7 = 00, while the latter provides 7 = 3. Hence, there is 
something wrong in their study. 

Their flaw lies in the two assumptions that were bor- 
rowed from those of Simon for the words problem in 
which the probability that the adding word is a new word 
is given /? while the probability that the adding word is 
the word that already exists i-times is given 1-/3(2^. 
In the word problem, there is only a single PA of word 
to the system and therefore, there is no difference in the 
sense of attachment of the new word and existing words 
since both are words. On the other hand, when we apply 
the method to a network model, we have to be careful be- 
cause in the network model, there are both attachments 
of nodes and links. We should not confuse with such at- 
tachments so that a node is exerted into the system at 
each time step and at the same time links are provided 
accordingly. We cannot separate both procedures using 
a single parameter f3 so that the attachment of a node is 
given by /3 and the attachment of links by 1 — /3. This 
is unacceptable. However, the mathematical ideas that 
they adopted is very interesting in its own right [29l|. 

In this paper, we would like to study the growing net- 
work models with PA with fractional power, using Si- 
mon's method [2^ in order to unify most of the previ- 
ously known methods and results studied in the litcra- 

Jp 111, Q E li6|. Ill liM [ii, & ill 

2T\ and to give new results. We would like 



to shed a new light on the problem. 

The organization of the paper is the following. In 
Sec. II, we review the growing network models with vari- 
ous types of PA. In Sec. Ill, we introduce the rate equa- 
tions for the connectivity distribution functions for the 
entire network and a particular node, respectively. In 
Sec. IV, we solve the rate equation for the connectivity 
distribution function for the entire network. In Sec.V, 
we solve the rate equation for the connectivity distribu- 
tion function for ith node. In Sec. VI, a conclusion is 
made. 
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II. THE EVOLVING NETWORK MODELS 

Random network models of Erdos and Renyi[l| have 
been standard models for complex network systems, until 
evolving network models have been found by Barabasi 
and Albert 0]. In the former, the size of the network N 
is fixed and the links are distributed with probability p 
such that each link has the mean connectivity (fc) = pN . 

In the latter, the network grows from a set of mo seed 
nodes, putting a newly added node from which m new 
links are distributed to the existing nodes. The way of 
adding a link depends on the PA of 



n, 



(1) 



This yields an evolution equation for the link number ki 
of the ith node: 



dk'i TTlhj'i ki 

—— = mlli = — jrr-, = — , 

dt Y.Uk^ 2V 

N-1 



(2) 



where we have a trivial relation ^ '2mt — m. 

Recently, Liu, Lai, and Ye[l^ have generalized the 
above PA to the following: 



n = ^ 

and the evolution equation to the following: 



dki 
~dt 



mil, 



mfcf 



k? 



(3) 



(4) 



where < a < 1, m is the average number of attached 
links and we have used the relation 



N 
3 = 1 



k" = II f 
r\jj — Pa'': 



(5) 



where 1 < /ic < 2m for < a < 1 and /^i = (fc) = 2m 
for a = 1. The solution of the above evolution equation 
has been given by them. 

This formalism is basically the same as that of 
Krapivsky, Redner and Leyvraz[ll| as discussed by Liu, 
Lai, and Ye[l9| in their appendix, since the former solved 
the time evolution of Eq.(4) while the latter solved the 
rate equation of the connectivity distribution function 
which will be discussed in the next section. Further- 
more, it is also almost equivalent to the formalism of 
Dorogovtsev, Mendes and Samukhinp^ and Bianconi 
and Barabasi flT'l. 

To be more specific to this point, the PA for the model 
of Dorogovtsev, Mendes and Samukhinp^ is defined as 
follows: 



A + h 



(6) 



where A is called the initial attractiveness of node i, 
which describes the ability of the ith node to receive links 
from other nodes, i.e., the probability for young nodes to 
get new links. On the other hand, the PA for the model 
of Bianconi and Barabasi [ITj is given by 



n. 



(7) 



where rji is called the fitness of node i, which describes 
the discrimination of different types of nodes. Although 
in the original Bianconi-Barabasi model the fitness 
is distributed at random with a distribution of p{r]) , we 
restrict ourselves to the case of the fixed fitness for each 
node. In this sense, we would like to call such fixed value 
of fitness the preference throughout this paper. Thus, we 
are able to recognize that the concepts of attractiveness 
and preference are essentially equivalent to each other, 
since one can expand the preference as rji — 1 + A/ki. 

In this way, we are led to define the most general form 
of the PA for a finite system of N nodes: 



n, 



(8) 



where a can be any real number of — oo < a < oo in 
general. 



III. THE RATE EQUATION 

A. Definition of tiie Rate Equation 

Let us first define that at each instant of time a node 
is added to the system such that the total number of 
nodes is given by N{t) — t. Let us denote by Nk{i, t) the 
distribution of the connectivity (= degree) k at the node 
i and define Nk{i,i) — 6k,m, which means that the node 
i starts with fci(O) = m when it is born[30j. 



Following both Krapivsky, Redner and Leyvraz ll| and 
Dorogovtsev, Mendes and Samukhinp^. we can define 
the rate (or master) equation for the time evolution of 
the connectivity distribution function as follows. Let us 
define the probability that the node i receives exactly I 
new links of the m injected as 



Pk,i 



7 ) (n,)'(i-n,^"-' 



(9) 



where 11^ is the PA given in the previous section, regard- 
ing ki as k. Then the rate equation is given by 



(10) 



1=0 



We would like to note here that Eq.(9) describes that 
at each time step, each node can accept I new links from 
a same single node, in general. However, in the most 
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of models the number of such receiving hnks is or 
1. Therefore, at this moment the above setting seems 
somewhat strange but the selection of the receiving link 
number can be realized from expansion of the PA [see 
Eq.(13)]. 



B. The Rate Equation for the Entire Network 

Let us define the connectivity distribution of the entire 
network by 

1 * 

Nk{t)^-J2Nk{i,t). (11) 

1=1 

Summing up Eq.(lO) over i from 1 to t, we obtain 

rn 

{t+l)Nk{t+l)-tNk{t+l,t+l)^"^tpk,iNk-iit). (12) 

1=0 

Now we expand the right hand side with respect to PA 
Ilfc as 

rhs = t{l-Uky^Nk{t)+mtUk-i{^-Ukr-^NkMt)+- ■ ■ 
= t{l - mUk)Nkit) + mtUk-iNk-iit) + ■■■ 



t - a^] Nk{t) + ^"'-'^"-'^" Nk-^it) 

Ma 



where we have defined 



It is equivalent to 
1 * 

i=l 



k=l 



(13) 



(14) 



(15) 



where 1 < fj,a < 2m{r]k) for < a < 1 and the appear- 
ance of the distribution function Uk in the second expres- 
sion is due to the definition that Uk — \ ^k,ki ■ If 
7]k — 1, then by definition /i^ — (fc"). The verification 
of Eq.(15) will be discussed in the next subsection[see 
Eq.(30)]. Here, we would like to note that the Ha in 
Eqs.(14) and (15) is different from the one in Eq.(5), 
since in the latter there is the preference factor in the 
expression. 

From Eqs.(12) and (13), we obtain 



{t + l)Nk{t+l)~t5k, 



t- 



m,rjkk" 

Ha 



Nkit) 



where higher terms are omitted. Exchanging some terms 
in both sides, we have 

{t + l)Nkit + l)~tNk{t) 



= —[-Tjkk''Nk{t)+f]k-i{k-irNk-i{t)]+tSk^m. (17) 

Ha 

This is the most general form for the rate equation for 
our purpose here. In these equations, Nk{t) are defined 
for k > m. 

Now, supposing that the time is as sufficiently large as 
t >> 1 such that t + 1 » t, the left hand side becomes 
t[Nk{t +1) - Nk{t)]. Then dividing both sides by t, we 
obtain 

Nk{t + 1) - Nk{t) 

Tfl 

= —[-7^kk''Nk{t)+rjk-i{k-irNk-i{t)]+6k^„r- (18) 

Ha^ 

This is exactly the same form of the evolution equation 
that was studied by Simon[l6l.[28j. apart from the coeffi- 
cient of the last term in the right hand side; in the Simon 
model, the last term is (35k fl, where (3 is the probability 
that the newly adding word is a new word while 1 — /? 
the probability that it is one of already existing words. 



C. The Rate Equation for the Connectivity 
Distribution Function 

We follow the same argument for derivation of the 
rate equation for the connectivity distribution function, 
Nk{i,t). Let us expand the right hand side of Eq.(lO) 
with respect to /. We obtain 



Nk{i,t+1) = pk,QNkii,t)+pksNk~i{i,t) + 



(19) 



Expanding the pk^ and pk,i with respect to the Ilfc in 
the same way, we obtain 



Nkit,t+1) 



1 



mrjkk" 
Hat 



Nk{i,t) 



mrik-i{k - 1)" 
Hat 



Nk-iii,t) + 



(20) 



This is the most general generalization for the rate 
equation of the connectivity distribution functions. It 
includes all the models studied in the previous literature 
U 



U 

2t 



12M 



IV. THE SOLUTION OF THE RATE 
EQUATIONS FOR THE ENTIRE NETWORK 

A. Simon's method 



^ m,k-.ik-ir ^^_^^^^^ 

Ha 



(16) 



Let us now consider the steady state solution of the 
rate equation for the entire network of Eq.(17). Let us 
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first note that if we expand the rate equation with respect 
to time t, then we are able to obtain the continuous rate 
equation such as those studied by Krapivsky, Redner and 
LeyvrazfTT'l and Dorogovtsev, Mendes and Samukhinp^. 
Therefore, we can use their methods here. However, as 
is mentioned before, the rate equation has exactly the 
same form that was first studied by Simon for the words 
problem [28|. Therefore, we can apply his mathematical 
method to this problem as well. 

Following Simonj28l]. for the connectivity distribution 
function for the steady state we can invoke the following 
relation: 



Nk{t + 1) _ tj^ 
Nk{t) ^ t ■ 



(21) 



This mean that the Nk{t) grows like Nk{t) oc t. Next, let 
us define 



Nkit+1) Nk{t) _^ 

= Pk- 



Nk-i{t + l) Nk-i{t) 



(22) 



Substituting the relation of Eq.(21) into the left hand 
side of Eq.(18) and the relation of Eq.(23) into the right 
hand side of Eq.(18), we obtain 



t+ 1 



m 

fiat 



7;fc_i(fc- 1)° 



Nkit)+6k^m. (23) 



Hence, we can define the steady state solution that is 
independent of t such as 



Nkit)=nlt. 



(24) 



For fc > m, we can derive the following for the steady 
state solution nl: 



vk-iik-iy 



(25) 



For k ^ m, since Nk{t) = for < < to, Eq.(18) yields 



V(t)^l_^!!i!^7V„(t). 



Therefore, we have 



-t = n^t. 



(26) 



(27) 



Thus, from Eq.(25) we can derive for the steady state: 

k 



n*k = Pknl_i n,;, J]^ 



Vk-l ■■■Vn 



r(fc) 



r(m) 



{^+Vkk'^)---{^+Vmm'-)'- 



(28) 



where k > m and in the last step we have used the gamma 
function r(fc) = (A: - 1)! and Eq.(27). 

We can also rewrite the above in the same form that 
was introduced by Krapivsky, Redner and Leyvraz[ll| as 
follows: 



vk-i{k-ir 



Vkk" 



Vkk" ^ +?7„+i(to+ 



??m+l(TO+l)" rymTO" 



^+r]kk°' ^+r?„+i(TO + l)" 77fcfc" 



rikk°' 



nr 



^ P{k). 



(29) 



j—7n 'i^^VjJ 



Substituting Eq.(29) into Eq.(15), the definition for fia^ 
we have the self-consistency condition for 



oo oo k ^ 



(30) 



k—m 



k—mj—m "^?7jj 



Then by expanding the right hand side with respect to 
fia and subtracting some terms in both sides, we obtain 
the simphfied relation: 



^ oo k ^ 

m E n 1 + -I 



(31) 



k=m+lj-m+l 



This equation corresponds to Eq. (6) of Krapivsky, Red- 
ner and Leyvraz[ll| and if we assume rjk ~ 1, it coincides 
with Eq.(A.3) of Liu, Lai and Ye[l9| where the relation 
between a and /Xq, is demonstrated. 

If we use the general notation Ak for mrjkk" , then we 
may rewrite Eq.(29) as 



P{k) ^ nl 



A 



^ 1 

k n TT^' 

j—m ^3 



(32) 



which is equivalent to the Eq.(8) in Krapivsky, Redner 
and Levvraz[ll| if we take m = 1. From the definition of 
the PA, it is obvious that A^ is nothing but the functional 
of PA. Here we obtain 



^k-l{k-\Y 

l^+Vkk- 



+ '7m+i(m + l)" 



OO ^ 



(33) 



k—7n-\-l j— m+1 
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Eq.(29)[or Eq.(30)] include many cases. If a = 1 and 
rjk = 1, which is the case of Barabasi and Albert Q and 
produces = — '^ni, then 



P{k) ^ nl = 



2m(m + 1) 
k{k + l){k + 2) 



oc k-^. 



(34) 



If a = 1 and r]k ^ 1, which is the case of Krapivsky, 
Rcdner and Levvraz|ll| and produces Aoo — mrjoc, then 



P{k) = n'^ (X k , 



(35) 



where 7 = 1 + fii/Aoo = 1 + ^ii/ mrjoc and /ii can be 
calculated by using Eq.(30). 

Thus, the heart of Simon's method is equivalent to 
that of Krapivsky, Redner and Leyvraz[llJ. It provides 
a different result from that of Bornholdt and Ebelp^. 
This is due to the fact that both methods adopt different 
assumptions to make the network systems, although both 
use Simon's method. 



B. The Method of Dorogovtsev, Mendes and 
Samukhin 



Let us next consider how to solve the Eq.(18), once 
again. This time we follow and generalize the method 
of Dorogovtsev, Mendes and Samukhin|,12i]. Recently, a 
similar approach has appeared [26j. 

Going back to Eq.(17) and considering the steady state 
solution we have 



Nk{t) 



-77fcA:"iVfe(t) + Tjk-i{k - l)"iVfc-i(t)] 



(36) 

Substituting Nk{t) = nut into the above, then we obtain 



(37) 



For the sake of simplicity for our purpose here, we take 
rjk = I- Then Eq.(37) becomes 



m 



Uk - (fc - l)"nfc_i = —5k,m- 
m 



Now, let us define a generating function: 

00 

$,„(z) = UkZ^. 



(38) 



(39) 



k—tn 



Multiplying to Eq.(38) and summing up over k from 
TO to 00, we obtain 



— + (1 - z) U — 
TO \ az 



TO 



(40) 



where we have used the definition of the fractional 
derivativefalIsS,!!!: 



When a = 1, /ii = 2to and therefore, Eq.(40) turns out 
to be 



2 + (1 -z)z 



dz 



(42) 



This is essentially equivalent to Eg. (7) in the paper of 
Dorogovtsev, Mendes and Samukhin [12|. 

Let us show this. Define as $(z) = X^fcLo ~ 
n„i + Um+iz + • • • , which means that k starts not from m 
but from 0. Therefore, we can rewrite ^m{z) as ^m{z) = 
z™$(z). From this, we find 

Substituting Eq.(43) into Eq.(42), we obtain the follow- 
ing equation for <I>(z): 



z(l - z)— $(z) + (to + 2 - TOz)$(z) = 2. (44) 
dz 

By differentiation for both sides of Eq. (44) with respect 
to z, we obtain 

,d2$(z) , . . , d<i>(z) 

Z{l-Z) , !, ^ + TO + 3 - TO + 2 Z ;^-TO$(z) = 0. 

dz 



dz"^ 



(45) 

The solution of the above equation can be obtained by 
the Gauss' hypergeometric function: 



2Fi(a, 6, c;x) 



Tic) 



E 



T{a)T{b)^^ T{c + n)n\ 



(46) 

which is the solution of the following ordinary differential 
equation: 

, , d^y(x) . , , , , dy(x) , , . 
x{l ~ x) f ; ^ + [c - (a + 6 + l)z] - ahy{x) = 0. 

dx'^ 



dx 



(47) 

Using the boundary condition that $(0) = const and 
$'(0) = for Eq.(lO), we get $(0) = 2/(to + 2). Using 
this and with replacement of a = 1, 6 = to, and c — to +3, 
<f>(z) can be obtained as 



ci>(z) = 



TO + 2 



2Fi(l, TO, TO + 3; z) 



E 

fe=0 



2to,(to + l)r(A: + to) k 

r(fc + TO + 3) ^ ■ 



(48) 



Therefore, comparing it with the definition of <i>(z) 
Z]fc=o nk^mZ^, we find 



P(k) = Uk+r, 



2TO(m + l)r(fc + to) 

r(/c + TO + 3) 



d_ 

dz 



$m(z) 



E^^ 



(41) 



2to(to + 1) 



(fc + ■m){k + TO + l){k + TO + 2) ' 



(49) 
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Since in this expression k starts from fc = 0, if we change 
it to start from k — m (i.e., we change the variable of k 
from fc + to fc), then we obtain 



P(fc) = Uk 



2m{m + 1) 
fc(fc + l)(fc + 2) 



k-\ 



(50) 



Hence, our result above is exactly equivalent to Eq.(34) 
and coincides with that of Dorogovtsev, Mendes and 
Samukhin[12i]. We can repeat the similar procedure for 
the integer a case. 

Thus, we can understand that both methods of 
Krapivsky, Redner and Levvraz[T]| and Dorogovtsev, 
Mendes and Samukhinp^ give the identical result when 
a = 1 and, therefore, they are essentially equivalent to 
each other. However, this method is not directly ap- 
plicable to the case of fractional a. In this case we 
need a more sophisticated method such as the fractional 
calculus^)\. [s^. 33]. This is a really challenging problem 
for physicists. Some hints are discussed in Appendix A. 



V. THE SOLUTION OF THE RATE EQUATION 
FOR THE DISTRIBUTION FUNCTION 

In this section we are going to solve the rate equa- 
tion of Eq.(20), unifying both methods of Krapivsky, 
Redner and Leyvraz[ll| and Dorogovtsev, Mendes and 
Samukhin[13|. 



A. The Discrete Rate Equation 

Let us first apply the method of Krapivsky, Redner 
and Leyvraz[Tll|. Now, we want to solve the discrete rate 
equation: 

7Vfc(z, t 1) = A^fe(z, i) - ^^^A^fc(z, t) 

fiat 

^ m,,-,{k-ir ^^_^^^^^^^ ^^^^ 

Hat 

Let us define the average connectivity k(i,t) by 

oo 

k{i,t)=Y,kNk{i,t). (52) 

k—ni 

Multiplying Eq.(20) by k and summing up over fc, and 
using this definition we obtain 



k{i,t+l) = k{i,t) + —rikk°'{i,t), (53) 



where r]kk'^{i,t) = J2T^rn'nkk'^^k{i,t)- If we set 7]k = 1, 
then it becomes 

jji 

k{i,t + l) ^ k{i,t) -\ k°'{i,t), (54) 

flat 



where k°'{i,t) = J^'kLm ^"^kiht). Hence, if we can fur- 
ther impose a = 1, then we obtain 



k{i,t+l) = k{i,t) + ^^k{i,t), 



(55) 



where we have used fii = 2m. 

Eq.(54) can be regarded as a discrete version of the 
continuous time equation such as Eq.(2) with the PA Hj 
of fractional exponent a : 



dki ^ mkf 

—— = mill = r^-^ 



fiat 



(56) 



In general, since k^ii^t) is not equal to fcj = [k{i,t)]°', 
the discretization of Eq.(56) is not equivalent to Eq.(54). 
However, for a = 1, Eq. (56) becomes identical to Eq.(2) 
and therefore, the discretization of Eq.(2) becomes equiv- 
alent to Eq.(55). Eq.(55) can be solved directly as fol- 
lows: From Eq.(55) we find 

k{i,t + i) = h + lj k{i,t) 



1 



1 



1 



1 



2t J y ' 2{t-ti) 

Taking logarithm of both sides, we obtain 



fc(z, ti 



(57) 



lnk{i,t) = X! ( ^ + F 



In k(i,ti 



1^1 



1 . ft 



-^-+lnfc(z,i,)«^ln(^-j +lnfc(^,^,), (58) 

where we have assumed that t is sufficiently large such 
that t ^ ti. Hence, we approximately obtain 



fc(i, i) = m ( — 
tj 



1/2 



(59) 



where we have used k{i,ti) = m. This exactly coincides 
with the result of time development of the seminal work 
of Barabasi and Albert's continuous modelQ. 

To go further to the case of fractional a, the discrete 
model is not so convenient since Eq.(53) is not simple 
equation of k{i,t). Thus, for this purpose, we have to 
use the continuous rate equation. 



B. The Fractional Differential Rate Equation 

Let us now take a continuous time limit, by which we 
obtain the following continuous differential equation: 

l^t—Nkii,t) = -rfkk"Nkii,t)+rik-i{k~irNk-i{i,t), 
m ot 



(60) 



7 



where Nk{i,t) is the connectivity distribution function 
of node i at time t. For our purpose here, we will set 
Vk — V = constant, for the sake of simplicity. 

Let us define the Z-transformation of the connectivity 
distribution function Nk{i,t) by 



(61) 



where ^{i,t;z) = YX=Q^k+m(i,t)z'' . Multiplying both 
sides of Eq. (60) by z^ and summing up over k and using 
the fractional derivative that was introduced in the previ- 
ous section, we obtain the following fractional differential 
equation: 

-p-^Ol'^rrM, t; z) + {l- z) (z—j $,„(i , t; z) = 0, (62) 

where fia = This corresponds to Eq.(40) in Sec.IV 
and also it corresponds to Eq.(14) in Krapivsky, Redner 
and Ley vraz ll| . 



C. The Solution of the Rate Equation with Linear 
Preferential Attachment 



When a is fractional we meet the same difficulty that 
has been discussed in the end of the previous section. 
Therefore, let us restrict ourselves to consider the sim- 
plest case of linear PA with a = 1. In this case, fii = 2m 
such that /?! = ^ and therefore, Eq.(62) becomes 

1 (9 d 

iriTT^^rnii, t; z) + {l- z)z— $„(i, t; z) = 0. (63) 
Pi at az 

If we take rj = 1 this corresponds to Eq.(42) and if > 1 
then it corresponds to Eq.(14) in Dorogovtsev, Mendes 
and Samukhin[l^. 

Let us solve Eq.(63). Using the standard method 
for solving linear partial differential equation, it yields 
the following differential equations for the characteristic 
curve: 



dz 



z{l-z) t 



dt d^r. 







(64) 



where in the denominator in the last equation means 
that d^„i — such that $,„(i,i;z) — ci = constant. 
Since by definition Nk{i,i) = Sk,m, we have the initial 
condition 



which then gives 



C2, 



(67) 



l-z., 

Substituting Eq.(67) into Eq.(65), we obtain the relation: 

C2 



ci 



1 + C2 



(68) 



which is the characteristic curve that we seek for Eq. (64) . 
Since the solution moves as a point along this curve from 
the initial time t = ti, the solution of Eq.(63) is given by 



^rni'i,t;z) = ci{i,t;z) = 



C2{i,t;z) 
1 -I- C2{i,t;z) 



(69) 



Now, solving Eq.(66) for C2 and substituting it into 
Eq.(69), we obtain 



$„(i,t; z) 



l-z{l-{U/tr 



(70) 



Adjusting with our definition of Eq.(61), we finally obtain 
the solution of Eq.(62) as 



$(i,t;z) 



/3i 



i-z\i-{u/tr 



(71) 



Now expanding the denominator of Eq.(71) using the for- 
mula: 



1 



l-z)^ 



r(7)/fc! 



(72) 



we obtain 



Hht;z) = Y. 



T{k + m) ft 



fc=0 



r(TO)fc! V t 



m(3i 



(73) 

Comparing this with Eq.(61), we end up with the follow- 
ing distribution function: 



r(fc + m) / ti 



r(m)fc! V t 



m(3i 



-4 



/3i 



, (74) 



where J^kLo -^k+m{i,t) = 1 is obviously satisfied. The 
Eq.(74) corresponds to Eq.(16) in Krapivsky, Redner and 
Ley vraz [Il[and Eq.(15) in Dorogovtsev, Mendes and 
Samukhin[13|. 



$TO(i,i;z) = ci = z" 



(65) 



where Zi means the initial value of z at time ti = i. 
Solving the first relation yields 



l-z U,: 



01 



(66) 



D. Time Development of the Average Connectivity 

Applying Eq.(74) together with the following formula: 



7Z 



(75) 
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the time development of the average connectivity of node 
i is calculated as 



fc=0 



Pi 



(76) 



which corresponds to Eq.(59) for the case of the dis- 
crete rate equation. If there is no effect of the pref- 
erence so that 77 = 1, then /?i = 1/2; In this case, 
Eq.(76) recovers the same time development of the semi- 
nal work of Barabasi and Albert 0] once again. And also 
it corresponds to Eq.(16) in Dorogovtsev, Mendes and 
Samukhinp^. respectively. 

Notify Eq.(35) for the scale- free exponent of 71 for a = 



1, where 71 = 1 



1 



If we set 77 = rjoo then 



from the definition of Pi = rj/2 we find the relationship 
between 71 and /3i: 



1 



71 



1 



(77) 



This recovers Eq.(17) in Dorogovtsev, Mendes and 
Samukhinp^. As pointed out by them, this relation is 
exact for any finite preference of 77 as long as a = 1. 
However, we do not know yet whether or not this rela- 
tion holds true even for the fractional a such as (3a = 
1/(7. - !)■ 



VI. CONCLUSIONS 

In conclusion, we have studied an analytical method 
of how to obtain the connectivity distribution functions 
for the various growing networks with PA of fractional 
power. First, in order to unify the growing network mod- 
els in a general point of view, we have discussed the vari- 
ous types of PA [See Eq. (8)] . Second, we have presented a 
general version of the rate equations for the growing net- 
works with PA of fractional power [See Eq.(16)]. Here 
we have presented the rate equations for the connectiv- 
ity distribution functions for the entire network and for 
an arbitrary node, respectively. Third, we have presented 
the way of solving the rate equations for connectivity dis- 
tribution functions for the entire network from the point 
of view of Simon's methodpsj. We have shown that our 
method unifies both methods of Krapivsky, Redner and 
Levvrazjllj and Dorogovtsev, Mendes and Samukhinfl^ 
into the same mathematical framework. We have intro- 
duced an idea of fractional calculus that the rate equa- 
tion for the connectivity distribution function with PA of 
fractional power becomes the fractional differential equa- 
tion. We have presented a scheme to solve the problem 
in the appendix, where we have introduced a new type 
of transidental function that we call the upsilon function, 
T(a, b, c; 2;) [See A. 13]. Fourth, we have presented the way 
of solving the rate equations for connectivity distribution 
functions for the ith node, using the generalize method. 

Thus, we have unified some of the previously known 
methods and results of connectivity distribution func- 



tions for such growing networks into a single analyti- 
cal theory using fractional calculus in the spirit of the 
seminal Simon's method of word problem[2^]. In this 
context, we have emphasized that in spirit, Simon's 
method is identical not only to the method of Bornholdt 
and Ebel[lB| but also to both methods of Krapivsky, 
Redner and Leyvraz ll| and Dorogovtsev, Mendes and 
Samukhin^l2j. However, the former provides a differ- 
ent scale-free exponent from the latter, since the defining 
conditions for the growing networks are distinct to each 
other. We have been able to establish a new method 
of fractional calculus to the problem of growing network 
models with PA of fractional power of a. This type of 
the fractional differential equation [See Eq.(40)] seems to 
be not yet known in this field of growing network mod- 
els, which is different from the one that is known as the 
fractional diffusion equations etc in the fractional calcu- 
lus [31I, [12, HI] . Therefore, we believe that our approach 
may shed a new light on solving such a novel type of 
fractional differential equations in its own right. 
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APPENDIX A: FRACTIONAL CALCULUS 

For the noninteger case of a, we need a different way 
of solving the rate equation. We have to invent a new 
method. This case seems to require the method of frac- 
tional calculus [31I, m, [3^ . 

Let us see this point. Go back to Eq.(40), once again. 



+ n - z) ( z— 
m \ dz 



<i>,„(z) = f^z^ 
m 



(40) 



This type of differential equation has been called the 
fractional differential equation when the degree of the 
derivative is fractional, which is called fractional deriva- 
tive. Recently, there have appeared many applications of 
the concept to various areas of physics [3ll. [33. Issj . 

From this context, our problem is one of such applica- 
tion to physics as well. How to solve the equation such as 
Eq.(40) with fractional or irrational a would be a math- 
ematical challenge to physicists. 

Using the properties of fractional differentials that 
(il)" = let us rewrite Eq.(40) as 



$,n(^) = ^ 



(Al) 
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which corresponds to Eq.(42). Using the relation of 
Eq.(43), Eq.(A.l) becomes 



Ma , (dm 
TO \azz 



d 

Z— h TO 

dz 



a-l 



<^{z) 



Mq 
TO ' 

{A.2) 

which corresponds to Eq.(44). Operating both sides of 
Eq.(A.2) by z-^, it turns out to be the following frac- 
tional differential equation: 

f^z^^z) + {l-z)(z^+ m] $(z) 
m dz \ dz I 



+ [-TO + (to - l)z] l^z— + toJ $(z) 0, (A. 3) 

which corresponds to Eq.(45). Comparing this with 
Eq.(47), we now easily see that if a = 1, then Eq.(A.3) 
turns out to be Eq.(48). Thus, the problem is mapped 
to how to solve the fractional differential equation such 
as Eq.(40) or Eq.(A.2) when a is fractional. 

Let us denote as z = e". Then, Eq.(40) and Eq.(A.2) 
can be written as 



TO \du 



TO 



and 



ra \ du 



TO 



(A.5) 



respectively. Similarly, Eq.(A.3) is transformed into the 
following: 



Ma d 
TO du 



Hu) + (1 



du 



a+l 



$(m) 



+ [-TO+(TO-l)e"] — +TO $(w) = 0, M.6) 
\du J 

Let us now define the Laplace transformation: 

/>OC 

= *(s) = / $(M)e-""dM. [A.l) 



The definition of fractional derivative requires the follow- 
ing operation: 



C 



h TO 

du ' 



= (s + to)"*(s). 



(A.8) 



The proof of these relations is straightforward. Using 
these, the Laplace transform of Eq.(A.6) becomes 



Ma 



s*(s) + (s + to)"+'«'(s) - (s + to - l)"+'*(s - 1) 

-to(s + to)"*(s) + (to - 1)(s + to - l)"$(s - 1) = 0, 

(All) 

which yields a recursion relation: 



* s) = -IT 7 - !)■ 

w (s-l-TO " ^ ' 



Please compare this with Eq.(25) in Sec. IV. Since 
we take r/^ = 1 for our purpose here, it is clear that 
Eq.(A.12) in the language of the (continuous) Laplace 
transformation corresponds to Eq.(25) in the language 
of the (discrete) Z-transformation, where the continuous 
variable s -I- to in Eq.(A.12) corresponds to the discrete 
variable k in Eq.(25). 

Looking at Eq.(A.12), if we are able to Laplace trans- 
form Vl/(s) into the original function $(z), then the prob- 
lem can be solved. However, it is not so simple since 
Eq.(A.12) is a recursion relation for ^'(s) and the vari- 
able s is a continuous one. Therefore, we cannot simply 
represent \E'(s) in terms of s such as in the case of 
Eq.(25). Rather, it should be regarded as a functional 
relation. For example, it is supposed to be like a func- 
tional relation such as r(z + 1) = zr{z), which defines 
the Gamma-function r(z), where if z = rt, integer, then 
r(n-f 1) = n\. 

From this understanding, we would like to regard the 
relation (A. 12) as a defining relation for a new (unknown) 
type of function. Let us define the following function 
Ta{a, b, c; s), upsilon function of s, such that 

T„(a, 6, c; s) = ■^^tt?-— b,c;s- 1), (A.13) 
a + o(s + cj" 

where a, b, c are real numbers. From this, as suggested 
by Aomotofs^l, we find that the upsilon function can be 
represented in terms of an infinite product: 



1 



1 °° 1 
t TT t 

f c)" 1 + f ^-r-V 



Ta{a,b,c;s) b{s + c)<^ 11 1 + ° i 

^ ' ' ' ^ ^ 71—1 o (s+c+n) 



(^.14) 



The convergence of the upsilon function can be domi- 
nated by the function of an infinite product in the de- 
nominator. If we denote it by 



Us;x) = 11 [1 



n=l 



{s + n)" 



(^.15) 



The Laplace transforms of other terms are given as fol- 
lows: 



then from Eq.(A.14) the upsilon function can be obtained 

as 



£[e"$(u)] =*(s-l). 



iA.9) 



e" h TO $(u) 

\ du ' 



^ (s + TO-l)"5'(s-l). (A.IO) 



Ta{a, b, c; s) — b{s + c)"- Ja{s + c; a/b). 



(A.16) 



The analytical properties and convergence of this func- 
tion /q(s; x) may be obtained by the standard argument 
of entire functions 13511. 
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At this moment, we do not know much about the 
mathematical properties of the function T^ia, b, c; s)[3^. 
However, if we define this function, then we can formally 
represent the Z-transformcd function in terms of the 
upsilon function as follows. Adjusting Eq.(A.12) with our 
definition of (A. 13), our 'I'(s) can be written in terms of 
Ta{a, b, c; s) as 

M,(s) = T„(— (A.17) 

Now, we invert ^'(s) by the inverse Laplace transforma- 
tion: 

= — / *(s)e''"ds 

= T-/ T4^,l,m;s e^"ds, (A18) 
where c is a constant. Replacing as 2; = e", we finally 



obtain the following formula for $(z): 

$(z) = — / T„ ( ^, 1, m; s) z'ds. M.19) 
27ri J-oo+ci ^ "1 J 

This is supposed to be the solution of the fractional dif- 
ferential equation of Eq.(A.3). 

Next, if we can expand ^(z) in terms of z such as 

00 

$(z)=^T„(^,l,m;fc)z^ (A20) 

fc=0 

where 1, to; fc) denote the expansion coefficients, 

then comparing the definition of the Z-transformation: 
<i>(z) = X^fclo '^k+mZ^, we obtain 

Tifc+m = Ta ( — , 1,to; fc) . (^-21) 

V TO / 

Hence, our problem is solved. 
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